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/(-R)-theories of gravity are reviewed in the context of the so called gravitational baryogenesis. 
The latter is a mechanism for generating the baryon asymmetry in the Universe, and relies on the 
coupling between the Ricci scalar curvature R and the baryon current. Gravity Lagrangians of the 
form C(R) ~ R n , where n differs from 1 (the case of the General Relativity) only for tiny deviations 
of a few percent, are consistent with the current bounds on the observed baryon asymmetry. 



The discovery of the accelerated expansion of the Universe [1] has motivated the developments of many models of 
gravity. These models are built up, typically, either in the framework of the conventional General Relativity or in the 
framework of its possible generalizations or modifications. In the last years, among the different approaches proposed 
to generalize Einstein's General Relativity, the so called /(-R)-theories of gravity received a growing attention. The 
reason relies on the fact that they allow to explain, via gravitational dynamics, the observed accelerating phase of the 
Universe, without invoking exotic matter as sources of dark matter. The gravity Lagrangian for these theories is a 
generic function of the Ricci scalar curvature R, and the corresponding action with the inclusion of matter reads 

S = ^2 J d i xy/^gf{R) + S m [g^,ip m ], (1) 

where k 2 = 8irG. It is a difficult ask to deal with higher order terms in the scalar curvature, thus the forms of f(R) 
that have been most widely studied in literature are f(R) ~ R + eR m and f(R) ~ R n [3,2,4-6]. A careful and relevant 
study about the possible form for /(i?)-Lagrangians has been performed by Olmo in [7]. There it has been shown 
that solar system experiments provide strong constraints on the possible dependence of f(R). 

-2A < f(R) < R- 2A + l 2 R 2 /2, (2) 

i.e. the gravity Lagrangian must be nearly linear in R, with non linear terms bounded by quadratic terms. Therefore 
the result (2) excludes non linear terms growing at low R [8]. In (2) A is the cosmological constant. 

The aim of this paper is to show that f(R) -theories of gravity provide a framework in which the gravitational 
baryogenesis may occur, and yields the observed baryon asymmetry in the Universe. The origin of the baryon number 
asymmetry is still an open problem of the particle physics and cosmology. Big-Bang Nucleosynthesis (BBN) [9] and 
measurements of CMB combined with the large structure of the Universe [10] indicate that matter in the Universe is 
dominant over antimatter; the order of magnitude of such an asymmetry is 

2*z»*<flio-u 
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where iig (ng) is the baryon (antibaryon) number density, and s the entropy of the Universe. As argued by Sakharov, 
the baryon asymmetry may be (dynamically) generated if the following conditions are satisfied [11]: 1) processes that 
violate baryon number; 2) C and CP violation; 3) out of the equilibrium. The last point is important because to have 
different numbers of baryons and anti-baryons is different, the reaction should freeze before particles and antiparticles 
achieve the thermodynamical equilibrium. However, as shown in [12], a dynamically violation of CPT may give rise 
to the baryon number asymmetry also in a regime of thermal equilibrium. 

Recently, within Supergravity theories [13], Davoudiasl et al. [15] have proposed a mechanism for generating the 
baryon number asymmetry during the expansion of the Universe by means of a dynamical breaking of CPT (and CP). 
In this approach the thermal equilibrium is preserved. The interaction responsible for CPT violation is given by a 
coupling between the derivative of the Ricci scalar curvature R and the baryon current [16] 

^ J d^x^—grd^R, (3) 

where M* is the cutoff scale characterizing the effective theory. For different approaches also based on gravitational 
baryogenesis, see [14]. 

If there exist interactions that violate the baryon number B in thermal equilibrium, then a net baryon asymmetry 
can be generated and gets frozen- in below the decoupling temperature Trj. From (3) it follows M^ 2 (d IJ ,R)J fi = 



1 



M~ 2 R(ng — ng). Therefore the effective chemical potential for baryons and antibaryons is ps = R/M 2 = —fig, 
and the net baryon number density at the equilibrium turns out to be (as T > ms, where ms is the baryon mass) 
rig = gbPsT 2 /6. gb ~ O(l) is the number of intrinsic degrees of freedom of baryons. The baryon number to entropy 
ratio is [15] 



Tig 15(? g R 
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where s — 27r 2 g* s T 3 /45, and g* s counts the total degrees of freedom for particles that contribute to the entropy 
of the Universe. g» s takes values very close to the total degrees of freedom of effective massless particles <?*, i.e. 
g* s ~ <7» ~ 106. ns/s is different from zero provided that the time derivative of the Ricci scalar is nonvanishing. In 
the context of General Relativity, the Ricci scalar and the trace T g of the energy-momentum tensor of matter {Tg 1 ') 
are related as follows 

R = -8nGT g = -8nG(l - 3w)p, 

where p is the matter density, p the pressure, w = p/p, and T g = Tg . R is zero in the radiation dominated epoch of 
the standard Friedman- Robertson- Walker (FRW) cosmology, characterized (in the limit of exact conformal invariance) 
by w = 1/3. This way no net baryon number asymmetry can be generated. However, as we shall see, a net baryon 
asymmetry may be generated during the radiation dominated era in /(i?)-theories of gravity (see [15,17] for other 
scenarios) 

The variation of the action (1) with respect to the metric yields the field equations 

f'R^ - | - V„V„/' + g^uf = k 2 T 9IXV , (5) 
where the prime stands for derivative with respect to R. The trace reads 

3D/' + f'R - 2/ = k 2 T 3 . (6) 



In the spatially flat FRW's metric 
Eqs. (5) and (6) become 



ds 2 = dt 2 - a 2 (t) [dx 2 +dy 2 +dz 2 }, (7) 



-^f'-{-+i-f"R = n 2 p, (8) 
a 2 a 



- + 2^ /' + J - - 2- f'R f"R 2 - f'R = k 2 p , (9) 
a a z J 2 a 

3f"'R 2 + 3 f'R + Q^f'R + f'R - 2/ = k 2 T 9 (10) 

where the dot denotes the derivative with respect to the cosmic time, p — T^ , pS* = T*j, and T g is the trace of the 
energy-momentum tensor, T g = p — 3p. Moreover, the Bianchi identities give a further condition on the conservation 
of the energy 

p + 2,- a ( P + p)=Q. (11) 

To derive a solution to Eqs. (8)-(ll), an explicit form for f(R) has to be specified. Following the models proposed in 
literature, we make the ansatz 

w) =(!)"> ( i2 ) 

where A is a constant with dimensions m p ~ 2 ^ 1 (mp ~ 1.22 10 19 GeV is the Planck mass). By using Eq. (12) and the 
ansatz a(t) ~ t a , Eqs. (8), (9) and (11) imply 

n = 2a. (13) 
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FIG. 1. Plot of the function g a vs a. The range in which the function g a is positive, hence T is defined, is 0.155 < a < 1/2. 

One can easily check that the trace equation (10) is fulfilled. 

The interaction (3) generates a net baryon asymmetry provided R ^ 0. As already noted, in standard cosmology 
R vanishes during the radiation era, deviations from the standard General Relativity prevent the Ricci curvature, as 
well as its first time derivative, to vanish. Directly from the definition of Ricci scalar curvature 

' a a 2 
R=-6[- + -z 
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it follows 



R 
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The net baryon asymmetry (4) turns out to be 



riB gb 45 a(l — 2a) 
~ ~ gZ^ 2 t 3 D T D M, 



(14) 



(15) 



where tv is the decoupling time. Equating the expression of p given by Eq. (8) to the usual expression of the energy 
density [19] 
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one obtains 



where 



1/2 



\4n 3 gJ 9a t a A<*/ 2 



(16) 



(17) 



2Q , a -l0a 2 + 8a-l 
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Fig. 1 shows g a vs a. The allowed range for the parameter a is 0.155 < a < 1/2 (the value a = 1/2 corresponds 
to Einstein's theory of gravity). Nevertheless, as we shall discuss below, BBN imposes stringent restrictions on its 
possible values 

Inserting (17) into (15), one has 



n B 



< 5.8 10- 3 # Q 



(18) 



where 



Hr, = 



y/a\l - 2a\ 



876,5- 



2(1 - 2a) 
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As pointed out in [15], a possible choice of the cutoff scale M* is M* = mpij\f%T; if Tu = Mj, where Mj <~ 2 10 16 GeV 
is the upper bound on the tensor mode fluctuation constraints in inflationary scale [18]. This choice is particular 
interesting because implies that tensor mode fluctuations should be observed in the next generation of experiments. 
For A <~ trip , using the constraint on the observed baryon asymmetry ub/s < 9 1CP 11 , Eq. (18) can be recast in 
the form 



A a <1.6, 



(19) 



where 



A Q = 10 12 # Q (2.75 l(T 3 ) 3/o 



(20) 



Eq. (19) is the main result of our paper: It relates the estimation of the baryon number asymmetry inferred from 
observations to the parameter a, hence to the exponent of the Ricci scalar curvature entering in the generic function 
f(R). This allows to fix the form of the /(i?)-Lagrangian. 

Let us now discuss /(i?)-theories of gravity in the context of BBN. According to BBN, the formation of light 
elements in the early time occurred when the temperature of the Universe was T < 100 MeV and the energy and the 
number density were dominated by relativistic particles. At this stage of the evolution of the Universe, the smattering 
of neutrons and protons does not contribute in a relevant way to the total energy density. All these particles are in 
thermal equilibrium owing to their rapid collisions. Besides, protons and neutrons are kept in thermal equilibrium by 
their interactions with leptons 



n 

( ' + n 
n 



p + e 

p + e~ +u e 



(21) 
(22) 
(23) 



The weak interaction rate A(T) is determined by means of the conversion rates of protons into neutrons, X pn (T), 
and the inverse ones X np (T). The rate X np (T) is expressed in terms of the sum of rates associated to the individual 
processes (21)-(23) [20,19] 



A 



np 



x 



+ A„+ e 



*p+v t 



■p+e +v e 



(24) 



1.293MeV, and 



Similarly for X pn (T). At enough high temperatures (T > Q and T 3> m e , where Q = m n — m 
m n p e stands for the mass of the neutron, proton and electron, respectively) the weak interaction rate is given by 
[19,20] 



A(T) = A„ P (T) + A P „(T)~ 



7tt(1 + 3g 2 ) 
60 



GlT 5 . 



(25) 



where Gf — 1.166 10~ 5 GeU~ 2 is the Fermi coupling constant and gA — 1-26 is the axial-vector coupling constant of 
the nucleon. To compute the freeze-out temperature Tf (the temperature at which baryons decouple from leptons), 
one has to equate A(T) given by (25) to the expansion rate of the Universe H = a/ a: A(T) ~ H. Using (17) with 

replaced by g[ BBN ^ = 10.75, it follows 



6.4 1 28 - 



.19 a 



1/4 



15 



GeV. 



(26) 



Coming back to Eqs. (19) and (20), the analysis of the gravitational baryogenesis indicates that the value a ~ 
0.15506 is compatible with the estimation (19). Nevertheless, such a result is inconsistent with BBN. As matter of 
fact, the plot of Tf vs a in Fig. 2 shows that values of a close to 0.155 have to be discarded because they lead to 
temperatures incompatible with the typical temperatures of BBN ((0.1 — 100)MeV). This is not the case for those 
values of a close to 0.5, as arises from Fig. 3. Thus we shall be interested only in the gravitational baryogenesys for 
a < 0.5. The results are reported in Fig. 4. As we can see, a narrow region for the parameter a, 0.46 < a < 0.49, 
is consistent with the order of magnitude given in (19). In particular, the maximum value of the function A a , 
^jnax) ^ i 59^ j s m excellent agreement with the estimate (19). Ai™ 02 ' 1 corresponds to a ~ 0.4845 so that n ~ 0.97. 
The gravity Lagrangian therefore reads C ~ i? 0,97 , i.e. the form of the function f(R) is nearly linear in R. We also 
note that in the limit a — > 1/2, the function A Q (hence H a ) vanishes and no net baryon asymmetry is generated. 



4 





A„ 



1.5 
1.25 
1 

0.75 
0.5 
0.25 




a 



0.35 0.4 0.45 0.5 
FIG. 4. Plot of A Q vs a. The dashed line represents the upper bound 1.6. The maximum value ~ 1.59, corresponding to 
a ~ 0.4845, agrees with the estimate (19). Moreover, in the limit a — > 1/2, A Q vanishes, so no net baryon asymmetry is 
generated (in agreement with General Relativity). 
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This is as expected because a = 1/2 implies that a(t) evolves as i 1 / 2 , i.e. the conformal symmetry is restored, and 
therefore the Universe expands according to General Relativity field equations. 

In conclusion, /(i?)-theories of gravity provide a natural arena in which the baryon asymmetry in the Universe may 
be generated through the mechanism of the gravitational baryogenesis. As we have shown, deviations from General 
Relativity, although very tiny (a few percent), allow to account for the observational data on the matter-antimatter 
asymmetry. These results rely on the exact solution of /(i?)-ficld equations (Eqs. (8) and (9)), which allows to 
determine (without a fine-tuning) the exponent n for the generic function f(R) ~ R n . The results are also consistent 
with BBN. 
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